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$H . Abstract. We discuss in this work the validity of the theoretical solution of the nonlinear Couette flow for a granular impurity 

Oh' obtained in a recent work [preprint arXiv:0802.0526|, in the range of large inelasticity and shear rate. We show there is a 

•^r ' good agreement between the theoretical solution and Monte Carlo simulation data, even under these extreme conditions. We 

I also discuss an extended theoretical solution that would work for large inelasticities in ranges of shear rate a not covered by 

>0 , our previous work (i.e., below the threshold value Oth for which uniform shear flow may be obtained) and compare also with 

simulation data. Preliminary results in the simulations give useful insight in order to obtain an exact and general solution of 

the nonlinear Couette flow (both for a > flth and a < Oth)- 
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j^ . Introduction. Granular materials are generically characterized by the inelasticity in the collisions between the par- 

ri ' tides they are composed of. Transport theories of fluidized granular materials have been the subject of a considerable 

I , amount of research work in the last few years (11121 ■ We may find two important reasons motivating this effort: 1) the 

'^ ■ transport of granular fluids has numerous and obvious industrial applications, and 2) the theories and experiments on 

G I granular media constitute, in several important cases, generalizations of existing theories on elastic fluids 1 1, 2]. Ad- 

ditionally, computational tools used for elastic fluids such as molecular dynamics (MD) |3] and Monte Carlo (DSMC) 
simulations |4] may also be used in the realm of inelastic collisions 111,|2D, and for this reason the study of granular 
media is an open and challenging field also for computational physicists interested in non-equilibrium systems [ 1 ] . We 
may illustrate the idea of generalization of previous theories with the help of a paradigmatic example: in the smooth 
hard sphere model, inelastic collisions are characterized by the coefficient of normal restitution a, whose value in the 
particular case of elastic collisions is unity L2] . For instance, the recent calculations of transport coefficients of granular 
gases for a monocomponent gas Q, binary dilute mixtures la], and the more realistic case of multicomponent dense 



^SJ ' mixtures [7] are based on the extension of the Chapman-Enskog method to inelastic gases, and this is done by taking 

into account inelasticity (i.e., the possibiHty of having a ^ 1) in the collisional integrals of the coiTesponding kinetic 
^\ • equation |l5|l ■ The multiplicity of behaviors of granular systems 1 1 ] emerges as a very important consequence of the 

QQ I fact that the study of granular gases is a generalization of that for elastic gases: for example, while for elastic gases the 

^^ . steady Couette flow is characterized by a parameter 7 > 0, for granular gases the three possibilities 7<0,7=0,7>0 

can occur (we will later describe the meaning of this parameter), which results in having up to 6 types of steady Couette 
flows in granular gases whereas in elastic gases only one is possible |8]. Of course this complexity in granular media 
arises also in their rheology, whose study is particularly interesting because nonlinearity is inherent to granular rapid 
flows ISi], except in the quasielastic limit |8]. We will analyze in this work the rheology of a granular impurity under 
Couette flow, focusing on conditions of large inelasticity and hydrodynamic gradients. We will use for this the exact 
solution of the nonlinear Couette flow of a granular impurity, for the case 7 > 0, recently developed by the authors 
PO] and we will show, with the help of DSMC data, that our theoretical solution works well even under these extreme 
conditions. Nevertheless, our theoretical solution only covers the cases a > Oth, where a± is the value for which the 
well known uniform shear flow may be obtained [SJ], while the real situation is not restricted to this possibility, as sim- 
ulation data we present now clearly show. For this reason, we will also explore in this work the possibility of extending 
the solution to the case 7 < 0, which occurs for a < Oth. 

Results and discussion. This work is motivated by a recent previous work fw\, where we obtained the exact 
solution of the nonlinear Couette flow of a granular impurity for a BGK-type kinetic model 1 11], suitably adapted to 
the granular binary mixture with the addition of a drag volume force 1 12]. We assumed in our previous work ][1P] that 
the impurity (species 1) has the same flow velocity as the excess component (species 2), i.e., ui = U2, and that the 
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FIGURE 1. (a) The curvature parameter j for the excess component as a function of the reduced shear rate squared a^. The 
vertical dashed line indicates the threshold value cq^ separating the regions with y < (to the left) and y > (to the right). The 
continuous line stands for the theoretical solution, and the points for DSMC data. The extended theoretical solution for y < works 
reasonably well in a finite interval of a to the left of a^^^. (b) The temperature ratio x vs. cP' for a granular impurity with a = 0.5, 
a unit size ratio, and mass ratios mi jm^ = 2 (continuous line, solid symbols) and m\ jm^ = 0.5 (dashed line, open symbols). Lines 
and symbols stand for theoretical and DSMC data, respectively. The behavior of % clearly indicates energy non-equipartition even 
without shear (a ^ 0). 

concentration ratio n\jn2 and the temperature ratio T\ /Ti are spatially uniform. In addition, the hydrodynamic profiles 
of the excess component \y\ are an extension of those characterizing the nonlinear Couette flow for the elastic gas 
lUin . Thus, the hydrodynamic profiles have the following form: 
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-7= const, n272 = const, 
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where the x axis coincides with the direction of the flow, the walls of the system are parallel to it and perpendicular to 
the y axis, 7 is the curvature parameter for the excess component temperature, and V2 (y) "^ n2 (y) v^zOO i^ an effective 
collision frequency |10]. The solution provides 7, the temperature ratio x = T1/T2, and the rheological properties 
of both species as functions of the reduced shear rate a and the mechanical parameters (size ratio, mass ratio, and 
coefficients of restitution). In the solution use is made of the characteristic function /^o.m(7, C2/V2) (where ^2 is the 
cooling rate of species 2) given by jlOl] 
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Note that Fqi„{y, C2/V2) is well-defined for 7 > only. As said before, the case 7=0 corresponds to a threshold value 
ath of the shear rate at which the symmetric Couette flow reduces to the uniform shear flow [10,^. On the other hand, 
one can construct an analytical continuation of the function /t)_m(y,z) for negative y values as 
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We proved the hypotheses ui = U2, «i /n2 = const, T\/T2^ const, and Eq. dU to be accurate by numerically solving 
(DSMC) our model kinetic equations, at least for values of the coefficient of normal restitution as low as a = 0.8 
and values of the reduced shear rate as high as a ^^ 1 1 10]. We present in this work results for considerably larger 
inelasticities (a = 0.5) and shear rates (up to a = 1.75). Additionally, we present preliminary results, for the excess 
component, of an approximate solution based on Eq. (|3]l of the nonlinear Couette flow in the 7 < region that we 
are currently working out. As we see in Figs. [T]and|2] the results are quite satisfactory: in Fig. [iffl) we can notice 
that the extended solution works well for 7(0) in a finite range of values of a < Oth, although the solution evidently 
fails below a critical value of a, and that the regular solution for a > a^^ (7 > 0) shows an excellent agreement with 
DSMC data. It is also to be noticed the good agreement for the temperature ratio x — T1/T2, that, as it is known and 
we can see in the figure, in granular mixtures is different from 1 even in the absence of shear (a — 0). In Fig. |2]we 
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FIGURE 2. 

coefficients vs 

coefficient M\ vs. a"^ for the same granular impurity. Lines stand for theory and points for simulation data. 



(a) Viscosity T]i (continuous line, solid symbols) and normal stress difference A'l (dashed line, open symbols) 
a^ for a granular impurity with a = 0.5, a unit size ratio, and a mass ratio m\/m2 = 2. (b) Normal stress difference 



present the results corresponding to the stress tensor transport coefficients of the impurity. All of them show a very 
good agreement with the numerical solution: the non-Newtonian shear viscosity rji, and the normal stress differences 
A^i = {Pi.xx^Pi.yy)/pi7Mi = {Pi.zz — Pi.yy)/pi, withpi —iiiTi, whosc magnitude shows the importance ofrhcologlcal 
effects in this system. We can also notice in Figs.[T]and|2]that, independently of the limitation of our theory to the case 
a > Oih, the numerical solution of the kinetic equation (simulation data) extends beyond that limit, and the functions 
represented do not show any discontinuity in their behavior (except perhaps Mi) at a = Oth. 

Summarizing, we have shown in this work that the theoretical solution of the nonlinear Couette flow obtained in 
a previous work lUOll can be safely used for large inelasticities (a = 0.5) and hydrodynamic gradients (a ^ 1.75). 
We also checked that the hypotheses underlying the theoretical solution, as well as its predictions for the rheological 
properties, are fulfilled for this large inelasticity. We have also discussed an extension of this theoretical solution in the 
region of shear rates below the threshold value a < ai^, i.e., negative value of the curvature parameter 7. These results 
are also useful for a future work, where we will assess to which degree our BGK-type kinetic model is able to describe 
the results from the 'true' Boltzmann kinetic equations. 
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